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Abstract 

We re-derive and discuss two low energy theorems which relate the 
potential of a static quark-antiquark pair with the total energy and action 
stored in the flux tube between the sources. In lattice QCD these relations 
are known as Michael's sum-rules, but we give an essential correction to 
one of them. Then we relate the low energy theorems to the virial theorem 
for a heavy quark-antiquark bound state. Finally we compare the results 
for the flux tube formation, which have been calculated in the model of the 
stochastic vacuum, with the low energy theorems and obtain consistency. 
From this we conclude that the model describes the non perturbative gluon 
dynamics of QCD at a renormalization scale, where the strong coupling 
constant is given by a s = 0.57. 
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1 Introduction 



As it is well known there is no analytical tool comparable to perturbation theory 
for treating QCD processes involving large distances. The most direct approach 
is a numerical evaluation of the lattice regularized version of QCD. Nevertheless 
it seems useful also to try analytical approaches, based on model assumptions, 
which may reveal some of the underlying structure of the theory. The success 
of the QCD sum-rules makes it plausible that the properties of the nontrivial 
QCD-vacuum play an essential role for the infra-red behavior of the theory. 
The model of the stochastic vacuum (MSV) makes the crucial assumption that 
the complicated infra-red behavior of non-Abelian gauge theories can be ade- 
quately approximated by a cluster expansion of the stochastic process describing 
the quantized theory U f2| H • It has been shown that such a model leads to 
linear confinement. In order to make more specific and quantitative predictions 
one has to make more radical assumptions, the most radical one being that 
all higher cumulants can be neglected as compared to the two point function. 
Then the stochastic process is a Gaussian one, i.e. all higher correlators can 
be reduced to products of two point functions by factorization. We refer to 
|| for the description of the model and its applications. In this note we 
present some consistency checks by comparing results of the model for heavy 
quark- antiquark systems with low energy theorems. 

Our paper is organized as follows: In section ||] we re-derive some low energy 
theorems, in section || we show how these theorems relate to the virial theorem 
for a heavy quark-antiquark bound state. In section || we present the results for 
the string formation calculated in the MSV and in section [| we compare these 
findings with the low energy theorems and discuss them. 

Except for section ^ we work with an Euclidean field theory (EFT), where 
the square of the electric field has the opposite sign of the same quantity in a 
Minkowski field theory (MFT), whereas the square of the magnetic field has 
the same sign in both theories: 

-^EFT — "^MFTi -"EFT — "MFT ■ 

Here and in the following a sum over the N% — 1 color components is understood. 



2 The Low Energy Theorems 

In this section we consider pure gluon QCD and static quark-antiquark sources. 
We follow the procedure of Novikov, Shifman, Vainshtein and Zakharov (NSVZ) 
H| and redefine the unrenormalized matrix valued gluon field tensor by multi- 
plying it with the bare coupling constant qq: 

F> = <7oF$ . (1) 
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Q 

The QCD-action Sqcd then depends on a s o = only in the form 

Sqcd = ^— I d 4 x Tr [F^(x)F^(x)] . (2) 



Let us consider the vacuum expectation value of a Wegner- Wilson loop ||] 0] : 

< W[L) >=< Tr Vexp[-i J A^dx^] > . (3) 

Here L denotes a closed loop and V path ordering. In a somewhat symbolic 
notation we write this functional integral: 

< W[L] > = —J VAeM-SqcD] Tr V J A^dx^] , 

M = J VAexp[-S Q cD] ■ (4) 

Differentiating log < W[L] > with respect to — g7r * - we obtain: 
2 31og[< W[L] >] 

87rQ sO 7T- = 

oa s0 

< ^ > < J d 4 x Tr [F^(x)F^(x)] Tr V exp[-i J^A^dx^} > 
-< |d 4 x Tr [F^(x)F^(x)] > . (5) 

If the loop L is a rectangle with "spatial" extension —R/2 < 23 < i?/2 and 
"temporal" extension —T/2 < 24 < T/2 the right-hand side (rhs) of eq.(|5|) can 
be simplified further in the limit of large T. Let the length a be of the order of 
the correlation length of the gluon field strengths. Then for \x^\ > T/2 + a the 
rhs of eq.([|) is zero and for \xi\ < T/2 — a independent of 24. This gives for 
eq.(|5|) in the limit of large T: 

87ra -° aiQg[ lTo [L]>] = T< J d3xTr >* • ( 6 ) 

Here and in the following 24 is zero and the expectation value < . >^ shall 
denote the EFT expectation value in the presence of a static quark-antiquark 
pair at distance R, i.e. static sources in the fundamental representation, where 
the expectation value in absence of the sources has been subtracted. 

We define the potential V(R) by: 

V(R) ^ lim - l0g< ^ [L]> (7) 
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and thus obtain the low energy theorem: 



- W* o _l_2 = _ < J d 3 x Xr [F^( x )F^(x)} > R . (8) 

In eq.(^) and eq.(||) self energies are understood to be subtracted. 

Up to now we have used the bare coupling go and bare fields. If we use the 
background gauge fixing || we have 

g = Z~ l g , A M = Z-V 2 Af , Z g Z l l 2 = 1 , A„ = g Af = gA^ , (9) 

where g and A M are renormalized quantities. The renormalized squared field 
strength tensor is given by: 

Tr [F^(x)F^(x)] = Z-\ Tr [F$(x)F$(x)] , (10) 

where Z F 2 was calculated || in background gauge with Landau gauge fixing to 
be 

Z F 2 = 1 + \go^~ log Z 
2 dg 

and thus we obtain 

P- = zz F , . (11) 
We use eq.(S)-eq.(ll) to express eq.(B) in terms of renormalized quantities: 



dV(R) 1 

— a 



. M g < / d3x Tr i F ^( x ) F ^( x )] >* • ( 12 ) 

In order to evaluate the left-hand side (lhs) of eq.(|i~2|) we use standard renor- 
malization group arguments: The only scales are the renormalization scale \x 
and R, thus the potential V = V(R, (j,,as) satisfies on dimensional grounds: 

{ R M ~ ^t) v{11, /U ' as) = ~ v{11, ^ as) ' (13) 

But for V as a physical quantity the total derivative with respect to /i must 
vanish: 
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where 0(a s ) = ^^a s (^). We get from eq.(|l3|) and eq.(p^): 

^ = -M ( )+ (15) 

and we finally obtain with eq.(|T2|) the low energy theorem 

{V(H) + «^^} = \^<j d ^ Tr [F M ,(x)F^(x)] >i? (16) 
or equivalently 

(V(U) + = \^<j^ {E{ X f + 5(*) 2 ) > R . (17) 



The eq.(|T7D is for the case of a linear potential just one of the many low energy 
theorems derived by NSVZ ||. But there renormalization was only discussed 
to leading order. 

We have also the usual relation between the potential V(R) and the energy 
density Qqq(x): 

V(R) = J d 3 x < e 00 (x) >R=\< J d 3 x (-E{xf + B{x) 2 ) > R . (18) 

Note that we are in a Euclidean field theory, hence the minus sign of the electric 
field, and that on the rhs of eq.(|i~7D and eq.(|i"8|) we have renormalized composite 
operators. 



In lattice QCD eq.([L7|) and eq.([l^) are known as Michael's sum-rules [10]. But 
in the derivation of the action sum-rule (eq.(|T7|)) scaling of V(R) with R was 
not taken properly into account and hence the second term on the lhs of eq.(|D 
is missing in [10]. 



3 The Low Energy Theorems For A Quark- 
Antiquark Bound State 

In this section we work in Minkowski field theory of gluons interacting with a 
dynamical heavy quark field q. The symmetrized energy momentum tensor is 
given by: 



kA w kA ' w kA ' 

( \ - -IV i F .F*\l 4- 



©kA = 2 Tr [F Kp F p x] + -g K \ Tr [F^F^] , 



eJ A = l -q[ lK T> x +i x T> K )q. (19) 
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Here is the covariant derivative. 

Consider a quark-antiquark bound state, which is an eigenstate of energy and 
momentum: 

\M(p)> = \qq(p)>, 
<M{p')\M{p)> = 2p°(27rf5 3 {p'-p). (20) 

We have from Lorentz invariance 

< M( P )\e s y u m -(o)\M(p) >= 2p, Pu . (21) 

In the rest frame of the bound state the total energy is given by p R = 2m q + E R , 
where the binding energy Eb is the sum of the mean values of the potential 
and kinetic energies: 

E B =< V > + < T > . (22) 

With eq.(|l|) and 6 sym ' = g^Q 8 ^ we obtain for \E B \ <C m q : 

<M(p R )\QZ m -(0)\M(p R )> = <M(p R )\&*^(0)\M(p R )> 

= 8m 2 q + 8m q E B ■ (23) 

The kinetic part of the energy momentum tensor of the quark-antiquark bound 
state is With the kinetic momenta of the quarks, (p q ) 2 = {pq) 2 = rn q , we 
have 

4 < M(p R )\Q%(0)\M(p R ) >=< **pZ + 2 ^f± > . (24) 

PR Pq Pq 



This gives for B Q = g^Q% and 9^ : 

-L <M(p R )\e^(0)\M( PR )> 
z Pr 

1 ... ,.„Q 



2m q - < T > , 



F < M(p R )\Q$ (0)\M(p R ) > = 2m q +<T> . (25) 
z Pr 

More details on the derivation of eq.(^) and eq.(|25|) will be given in a future 
publication. Now we evaluate the gluonic part of the energy momentum tensor 
by subtracting the kinetic part (eq.(|25|)) from the total tensor (eq. (f23"[) ) . We 
obtain 

V <M{p R )\& G (0)\M(p R ) > = <V>+2<T>, 



2p 



2/> 



R 

y <M(p R )\eg(0)\M(p R )> = <V> (26) 
'r 
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with = g^&nu- I n order to relate the kinetic energy to the potential one 
we use the virial theorem: 

1 „ 9V 

<T>= 2 <R BR > ■ 



Eq.(p6|) then gives 

^<M(p R )\& G (0)\M(p R )> = <V > + <R^> , 
V <M( PR )\eg(0)\M( PR )> = <V> . (27) 



Using the trace anomaly of the energy momentum tensor [11] 



6 G = Tr [F^F^] 
we finally obtain: 

<V> + <R jR > = W R <M{PR) ^s TT [ F ^(°) F ^(°)]I M (^)>' (28) 
<V> = < M(p R )\eg(0)\M(j> R ) > . (29) 



It is straightforward to show that here we have again the relations (16) and 
(|1S|), but in eq.(|2"8|) and eq.(p9|) we have mean values as obtained by integration 
with the absolute square of the wave function of the quark-antiquark bound 



state. The integration over the space points in eq.fllq) and eq.(18) corresponds 
to the translational invariance of the quark-antiquark bound state in eq.(^) 
and eq.pJ). 



4 The Flux Tube In The MSV 



In this section we present the main results from a calculation of the flux tube 
between a static quark-antiquark pair in the MSV. For all technical details we 



refer to our previous publication [12|. 



The expectation values of the squared elements of the gluon field strength tensor 
in presence of the quark-antiquark pair minus the vacuum values are given by: 

A ^ n _ 4 < W[L] P a p(x) >-< W[L] >< P aP {x) > 

a P {X > = W P < W[L] > • (30) 
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Here P a /3(x) is a Wegner- Wilson loop over a plaquette in a, /3-direction centered 
at x with side length Rp. For Rp — ► we have: 



lim AF„g (x) = -2 < <^ Tr [F Q/3 (a;)F a/3 (x)] >r (31) 
iip — >0 

where there is no summation over a and (3. 



It was shown in [12] that by symmetry arguments < g 2 B 2 >r= 0. The square 
of the electric field perpendicular to the loop L (< g 2 E 2 _ >) is also practically 
not influenced but only the squared electric field parallel to L (< g 2 E 2 >) is 
affected by the static sources. In fig.((l]) we display — < g 2 E 2 >r as a function 
of the perpendicular distance x± from the loop L and of the position along 
the quark-antiquark axis X3. The lengths are given in units of the correlation 
length a of the gluon field strengths, which is fixed to a = 0.35 fm in the MSV 



12 . 



R = 4a, a =0.35fm R = 9a, a =0.35fm 

Figure 1: Difference of the squared electric field parallel to the quark-antiquark 
(23) axis (— < g 2 E 2 (23, x±) >r) in for different quark separations R. x± 
is the distance transverse to the 23-axis and the dots denote the quark positions. 



As can be seen directly from eq.(|17]) and eq.fllq) the squared electric and mag- 
netic field strengths are separately not renormalization group invariant, thus 
statements like < g 2 B 2 >p= or < g 2 E 2 _ >r~ are scale dependent. In |l2| 



we have used eq.(18) in order to determine at which a s the MSV is supposed to 
work. We have calculated the total energy stored in the flux tube and compared 
it to the potential V(R), which can be determined directly in the MSV Jjjpl. 



In fig. (S) we present this comparison using our fitted value a s = 0.57 [12]. 



Figure 2: The total energy stored in the field (dots) calculated with the en- 
ergy sum-rule ( |i~8| ) as compared with the potential of a q-q-pair obtained by 
evaluation of the Wegner- Wilson loop in the MSV [|IJ [Q] (solid line) 



5 Discussion 



We first consider the region where the potential is linear. Then eq.(|i~7|) reads: 

l -^-< Jd 3 x (E(x) 2 + B(x) 2 ) > R = 2 V(R) . (32) 

Since in the model of the stochastic vacuum the static sources do not modify 
the color magnetic field, i.e. < B 2 >r= , we have consistency between eq.(|i"8|) 
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and eq.(|32|) only if 



P/a a = -2 . (33) 



From the energy sum-rule (eq.fllq)) we have already obtained a s = 0.57, thus 
consistency of the MSV requires that eq.([33|) should be satisfied for this value 
of a.,. 



The /3-function can be expanded in a power series in a s : 
Po_ 2 _ P}_ 3 _ Pi 

■ 27r a s g 47r 3 



P( a s) = ~^Z a s ~ T~o< " TTT-^^s + ■■■ ■ (34) 



The n-loop expansion gives the series up to a" • The first two terms of the 
power series are scheme independent, the numerical value of the term propor- 
tional to Ug has been evaluated [jDJ in MS. For pure gauge SU(3) theory we 
have @: 

0o = U, A = 51, /3 2I ^ = 2857. (35) 
Numerical investigations of the running coupling constant in lattice gauge 



Monte Carlo simulations [15| indicate that the power series fl34| ) for j3{a s ) trun- 
cated at a* can be used up to rather large values of a s in their scheme if the 
coefficient Pi is fitted numerically with the result pz ~ 6270. 

From fig.© we see that our condition (|33|) yields a s = 1.14 using /3 on the 
one loop level, a s = 0.74 on the two loop level and a s = 0.64 on the three 
loop level in MS. This value is already very close to the value of a s = 0.57 
determined from the flux tube calculation using eq.(|i~8|). If we use for (3 the 
series (|34|) truncated at order and impose condition ( |33| ) for a s = 0.57 
and determine the coefficient pi accordingly we find 02 = 6250 (see fig.©), 
i.e. about twice the MS value. Curiously this agrees almost exactly with the 



numerically determined value of [15| quoted above. 



Figure 3: as a function of a s for different /5-functions 



We see from fig.(|2|) that for distances R/a < 1 (a = 0.35 fm) the MSV does no 
longer yield a linearly rising potential and hence we have 

V(R) + R^^^2V(R), R/a<l. 



Therefore the conditions ( |l8| ) and ( |32| ) cannot be fulfilled simultaneously for 
< B 2 >r= 0. This does by no means speak against the model, since in the 
derivation of the results in |l^] it was stressed that the spatial extension of the 
loop had to be large as compared to the correlation length of the field strengths 
in order to justify the evaluation of E 2 and B 2 . 
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6 Conclusions 



We have re-derived and discussed low energy theorems and sum-rules relating 
densities of squared field strengths to the potential of a static quark-antiquark 
pair. We found that for a dynamic quark-antiquark bound state these sum- 
rules are in a sense equivalent to the dilatation trace anomaly and the virial 
theorem. We have applied the sum-rules to quantities obtained in the model 
of the stochastic vacuum. This model was introduced as a phenomenological 
approximation to the complicated measure for functional integration over the 
slowly varying gluon field strengths. No use of the equations of motion has 
been made explicitly. The factorization hypothesis (Gaussian measure) enters 
crucially in the evaluation of the squared color fields. In view of the drastic 
approximations made it is gratifying that relation (|33|), which is a highly non- 
trivial consequence of the dynamics of the system, is so well fulfilled. 

The scale /i at which the MSV works was found to be the one where a s (fj,) = 0.57 



in [12] by using the energy sum-rule. With the help of the action sum-rule we 
now checked the consistency of the MSV at this scale and obtained an effective 
/^-function. 



As can be seen from eq.(|i~7|) and eq.(|i~8|) the squared field strengths B 2 and 
E 1 depend on the renormalization scale \i and hence on a. s {p). We may use 



relations ( |i8|) and (32) in order to predict the ratio of 



Q _ Jd 3 x<B(x) 2 > R _2 + P/a s 



fd 3 x<E{x) 2 > R 2-P/a a 

as a function of the strong coupling a s . The result is plotted in fig.@ (we have 
used there our numerically determined value of fa). This ratio can be checked 
in lattice gauge calculations. 

Figure 4: The ratio Q as a function of a s 

In this paper we have worked in pure gauge theory and in a theory with gluons 
and one heavy quark field. We leave the investigation of similar sum-rules in a 
theory including light quarks for future work. 

Acknowledgments The authors thank E. Laermann, C. Michael, H.J. 
Rothe and I. O. Stamatescu for discussions and especially M. Jamin for valuable 
comments and suggestions. 



References 

[1] H.G.Dosch, Phys.Lett.B190 (1987) 177 

[2] H.G.Dosch, Y.A.Simonov, Phys.Lett.B205 (1988) 339 



9 



[3] H.G.Dosch, Prog.Part.Nud.Phys.33 (1994) 121 

[4] H.G.Dosch, E.Ferreira, A.Kramer, Phys.Rev.B50 (1994) 1992 

[5] V.A.Novikov, M.A.Shifman, A.I.Vainshtein, V.I.Zakharov, 
Nucl.Phys.B191 (1981) 301 

A.I.Vainshtein, V.I.Zakharov, V.A.Novikov, M.A.Shifman, 
Sov.J.Part.Nucl.13 (1982) 224 

[6] F.J.Wegner, J.Math.Phys.12 (1971) 2259 

[7] K.G.Wilson, Phys.Rev.B10 (1974) 2445 

[8] B.S.de Witt, Phys.Rev.162 (1967) 1195,1239 

J. Honerkamp, Nucl.Phys.B36 (1971) 130, B48 (1972) 269 
L.F.Abbot, Nucl.Phys.B185 (1981) 189 

[9] B.Grinstein, L.Randall, Phys.Lett.B217 (1989) 335 

[10] C.Michael, Nucl.Phys.B280 (1987) 13 

[11] M.S.Chanowitz, J.Ellis, Phys.Lett.B40 (1972) 397, 
Phys.Rev.B7 (1973) 2490 
R.J.Crewther, Phys. Rev. Lett. 28 (1972) 1421 
J.Collins, A.Duncan, S.Joglekar, Phys.Rev.B16 (1977) 438 

[12] M.Rueter, H.G.Dosch, Zeitschrift f. Physik C in press, [hep-ph 9409"362] 

[13] O.V.Tarasov, A.A.Vladimirov, A.Yu.Zharkov, Phys.Lett.B93 (1980) 429 

[14] Review of Particle Properties, Phys.Rev.B50 (1994) 1297 

[15] M.Liischer, R.Sommer, P.Weisz, U.Wolff, Nucl.Phys.B413 (1994) 481 



10 



4 4 




R = 4a, a =0.35fm i? = 9a, a =0.35fm 



Figure 1: Difference of the squared electric field parallel to the quark-antiquark 
(2:3) axis (— < g 2 E^(xs,xj_) >r) in for different quark separations R. xj_ 
is the distance transverse to the £3-axis and the dots denote the quark positions. 
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Figure 2: The total energy stored in the field (dots) calculated with the en- 
ergy sum-rule (18) as compared with the potential of a q-q-pair obtained by 
evaluation of the Wegner- Wilson loop in the MSV [1][2] (solid line) 
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